Quantum control of spin-correlations in ultracold lattice gases 
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We demonstrate that it is possible to prepare a lattice gas of ultracold atoms with a desired non- 
classical spin-correlation function using atom-light interaction of the kind routinely employed in 
quantum spin polarization spectroscopy. Our method is based on quantum non-demolition (QND) 
measurement and feedback, and allows in particular to create on demand exponentially or alge- 
braically decaying correlations, as well as a certain degree of multi-partite entanglement. 
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Ultra-cold atomic gases trapped in optical lattices offer 
an unprecedented playground for studying the quantum 
phases of many-body systems [1 j. In particular, quan- 
tum states of ultra-cold lattice gases with spin degrees 
of freedom may be used to simulate quantum magnetism 
and to investigate physics relevant for our understand- 
ing of high-T c superconductivity [2]. While enormous 
progress has been made towards engineering such sys- 
tems, achieving the regime of high-T c superconductivity 
remains experimentally extremely challenging because of 
the low temperatures required pQ. 

In this context, quantum spin polarization spec- 
troscopy (SPS) [3] has emerged as a promising technique 
for detecting quantum phases in lattice gases via the co- 
herent mapping of spin-correlations onto scattered light 
in a quantum non-demolition (QND) measurement. In 
particular, spatially-resolved SPS that employs standing- 
wave laser configurations [4] allows direct probing of mag- 
netic structure factors and order parameters [5]. In this 
manuscript, we propose inverting the SPS scheme in or- 
der to prepare a lattice gas with a desired non-classical 
spin correlation function. Motivated by the experimen- 
tal demonstration of spin-squeezing via QND measure- 
ments [6 ], and by the recent extension of these ideas to 
unpolarized ensembles [7 , we demonstrate that a simple 
modification of the experimental scheme of Ref. [4 (illus- 
trated in Fig. [I]) allows for the on-demand preparation of 
lattice gases with arbitrary spin-correlation functions. 

Atom-light interaction. — We consider the interac- 
tion of atoms trapped in a one-dimensional optical lat- 
tice potential with a set of standing-wave pulses of near- 
resonant light with wave-numbers k p . The atoms are 

described by collective variables J a ^ = Y^LiJ^al •> where 
the index n runs over the n a atoms at lattice site i and 
a = x,y,z labels the components of the atomic spin op- 
erators with length j. With n s lattice sites, the total 
number of atoms is Na = ^s^a- The photons are de- 
scribed by collective Stokes operators S a with a = 1, 2, 3, 

defined as S a = -(a+,a^_)a a (a+,a-) T , where the a a are 
the Pauli matrices, and a± are annihilation operators for 



the spatial and temporal mode of the pulse with circular 
plus/minus polarization. The atom-light interaction for 
a single pulse is then described by the effective Hamilto- 
nian [8] 



Hp — Qp ^ ^ c i(kp)J z ,iSs 



(i) 



where Ci(k p ) = (1 + cos(2fc p 7*i))/2 describes the standing- 
wave intensity profile, and the coupling constants ft p de- 
pend on the probe detuning and intensity. Eq. (fil de- 
scribes a QND measurement that induces spin-squeezing 
of the J z component of the collective atomic mode 
J a (k) = J27=i exp(z/cr i )/ v /n^ with k = ±2k p . For 
multi-level alkali atoms, this effective Hamiltonian can 
be synthesized using multicolor or dynamical-decoupling 
probing techniques [6j [9] . 

Measurement & feedback — We model the inter- 
action using methods developed for treating the 
Gaussian dynamics of collective-variable systems pj 
HQ] , with the assumption that n a ^> 1. The 
full system is then described by the operators 
R m (k) = {J x (k), J y (k), J z (k), Si, S2, Ss} and the covari- 
ances f mn (fci, fc 2 ) = (Rm(ki)Rn(k 2 )+Rn(k 2 )R m (k 1 ))/2- 
(R m (ki)} (R n (k 2 )) . The dynamical equations for the 
covariances can be derived from the Heisenberg equa- 
tion of motion for the operators, where, in the small- 
angle regime, an operator changes as R m (k)( out ^ = 
#m(&) (in) - ir[R m (k)^ n \H p ], where r is the pulse dura- 
tion. We assume that the atomic and light variables are 
initially uncorrelated, and that the atomic covariances 
T^g (ki,k 2 ) = Va 7^ ft. For simplicity we also assume 
a uniform atomic filling factor. 

For an input Si -polarized pulse, the only covariances 
that change due to the pulse are 
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FIG. 1. Proposed experimental set-up. Atoms trapped 
in an optical lattice (black) are probed with a far-detuned, lin- 
early polarized standing- wave light pulse with wavevector k p . 
Afterwards, the probe beam is out coupled to a homodyne de- 
tector, where S2 is recorded. A projection-noise limited mea- 
surement induces spin squeezing, introducing quantum corre- 
lations among the atoms in spatial mode k = 2k p . Feedback 
is applied via optical pumping to set (J a (k)) = 0. Succes- 
sive spin components J a can then be separately squeezed by 
coherently rotating the atomic spin between measurements. 
The procedure is repeated for a set of wavevectors k p , with 
interaction strengths Q p weighted by the corresponding am- 
plitude of the cosine Fourier transform of the desired spatial 
correlation signature. 



tQ p Si ^n s j/S 
f zz (2k p ,k)/2 + 



where the coupling strength C p - 
and we define T zz (k) = F zz (0,k) 
T zz (— 2/cp, fc)/2. 

Detection of 62 then transfers the correlations de- 
scribed in Eqs. pi to the atoms. This can be modeled as 

a projection iw= r(° ut ) - r(° ut )(n 2 r( out )n2) MP r( out ), 

where MP indicates the Moore-Penrose pseudoinverse 
and II2 = diag(0, 0, 0, 0, 1,0) [7]. After the measurement, 
the atomic covariances are 



f,(out) / 
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Eqns. Q and Q imply that if the covariance ma- 
trix T a/ 3 is initially diagonal, then the only atomic co- 
variances changed by the interaction-measurement pro- 
cess are T zz {k\, fc 2 ), i.e. the measurement induces spin 
squeezing of the J z (2k p ) mode. Furthermore, the process 



is highly symmetric, preserving T aa (k, k') = T aa (k, —k') 
for k^k' and T aa (k, k') = r aa (fc', k) Vfc, k' . 

The orthogonal spin components J a (2k p ) can be suc- 
cessively measured by coherently rotating the atomic 
spin between measurements. To allow the measurement- 
induced squeezing to be repeated for each spin compo- 
nent, we require (J a (2k p )) =0, which allows us to avoid 
measurement induced back-action due to the Heisen- 
berg uncertainty relation (A J a (/i:i)) 2 (AJ^(fc2)) 2 > 
1^1(^1 + fe))| 2 - To obtain (J a (2k p )) = 0, the 
result of the measurement of 62 can be used as the input 
to an optical pumping feedback process: a weak pulse 
of near-resonant light at wavevector k p with an intensity 
proportional to slf^ will set (J z (2k p ) + J z (—2k p )) = 
0, and a second pulse with a half-period phase shift 
sets (J z (2k p ) - J z (-2k p )) = 0, so that (J z (2k p )} = 
(J z (—2k p )) = 0. As shown in Ref. [7], this feedback 
introduces spin noise oc TV 1 / 4 , which is negligible in the 
thermodynamic limit [TTj . 

Strategy — We now motivate a strategy that exploits 
the spatial dependence induced by Ci(k p ) in Eq. ([!]) to 
systematically manipulate the real-space spin-spin cor- 
relation function 

r a a(ri,r 2 ) = (Ja(ri)J a (r 2 )} - (J a (ri)) (J a (r 2 )) (5) 
= — exp(zfciri)exp(zfc 2 r2)f aa (fci,fc 2 ) 

ki,k 2 =l 

where we label the /c-vectors in the first Brillouin zone 
(BZ) from 1 to n s . Eqs. ^ and Q imply that the co- 
variances of a given collective mode J a (k) are only al- 
tered by a pulse with k = 2k p [12]. This suggests that 
we can manipulate r aa (n,r 2 ) with a sequence of pulses 
with wavevectors 2k p that cover the first BZ. Now we 
show that this can be done by an appropriate choice of 
coupling constants C p . 

We assume that we start with a completely mixed ini- 
tial state with covariances r^(ri,r 2 ) = T^S ri)r2 where 
r(°) = n a j(j + l)/3. Under the approximation that the 
covariances T aC t(ki, k 2 ) with k 2 7^ ±fci,0 remain small, 
Eq. ([5| becomes 



T zz (r u r 2 ) « — f^(0,0) + — 5^[cos(2fepri)+cos(2fepr 2 )]f^(2fcp,0) 



2 r 

+ ~Y1 \ cos ( 2k p( r i + r 2 ))f zz {2k p , 2k p ) + cos{2k p {n - r 2 ))f zz (2k p , -2k p ) 



(6) 



The spatial dependence is strongly dominated by the define the scaled coupling constants 
T zz (2k p , —2k p ) term, which, after the pulse with k p , 2 

changes as r z ° z ut \2k p , -2k p ) = r(°>(l - / p /4), where we / = T C ' (7) 
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FIG. 2. Numerical results for (a) exponential decay, and (b) algebraic decay, (i) First column: Real-space spin 
correlation functions. The lines from darker to brighter shades of blue (and in the main panels from top to bottom) are for 
optical depth d — oo, 300, 99 and 33. The d — oo data is plotted with open circles to allow comparison with the fitted curves 
(green dashed lines). The insets in (a.i) and in (b.i) are log- linear and log-log plots, respectively, where for clarity we subtract 
C(5r — > oo), extracted from a fit. In (a.i), without decoherence, the decay follows an exponential fit. Deviations at large 
distance become stronger with increasing decoherence (decreasing d), and C(5r — »• oo) cannot be reliably determined, yielding 
deviations from straight lines. In (b.i), the curves are straight lines for all values of d, and algebraic fits are very accurate, 
(ii) Second column: The covariances T(k, —k) closely follow the cosine Fourier-transform of the desired correlation signature 
Cdes(^), even f° r small d. Deviations occur only primarily at small k where the pulse strength C p is high, (iii) Third column: 
If the entanglement witness W(5r, <j)) is negative, the system diplays non-classical correlations. In (a.iii) and (b.iii), the red 
lines are fits to W(5r, 0) with exponential and algebraic decay respectively. Note that the orange contour line in (a.iii) marks 
W (5r, 0) 0. 



Assuming that these covariances are not changed by 
subsequent pulses [13 j, Eq. ([6| becomes 
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f**(0,0) 
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The spatial dependence is given by the final term, which 
is the cosine Fourier transform of the f p . This suggests 
the following strategy for manipulating the spin-spin cor- 
relations r zz (n,r 2 ): Let T zz (r u ri + Sr) = C des (5r) 
be the desired output correlation signature. To deter- 
mine the coupling strength C p that should be used for 
each wavevector k p in order to create Cdes(Sr), we ap- 
proximate f p by the inverse cosine Fourier transform of 
— 4Cdes(^ r )/r(°). Further, in Eq. ([7]), we replace the co- 
variances r zz (ki, k 2 ) in the expression for with the 
completely mixed values T^ ). Both approximations are 
valid for realistic experimental parameters. Now, we can 
solve Eq. Q for C p , 
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22 
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(9) 



where g p = § for k p = and g v — \ otherwise. 

The coupling strengths C p can be adjusted experimen- 
tally by choosing detuning A, intensity, and duration of 
the pulse appropriately. In fact, C p = ^ 'iV A^Vph^T / 'AA, 
where a is the on-resonance cross section for the probe 
transition, 7 the spontaneous decay rate, and A the cross 
section of the atomic ensemble illuminated by the probe. 
With a finite on-resonance optical depth d, C p is re- 
lated to the probability of spontaneous emission rj p via 
C p = y/drjp, giving a trade-off between coupling strength 
and decoherence. Decoherence due to spontaneous exci- 
tation by the probe pulse is included in the model, fol- 
lowing Refs. 13 [lOj [14] , by updating the atomic covari- 
ances according to f ^(fci, £2) = (1 — 2r/ p )t aa (ki, k 2 ) + 
2 Vp T^S ku _ k2 . 

Entanglement witness. — A special kind of correla- 
tion, entanglement, is particularly important in the con- 
text of quantum information processing and many-body 
systems [15]. To show that our proposal can create mul- 
tipartite entanglement, we derive an entanglement wit- 
ness for the multimode spatial correlations induced by 
the procedure described above. Generalizing the strat- 
egy of Refs. [16], we use the witness W = S/n a — 1, such 



4 



that W < implies entanglement, where we define 



S = E 5 « = E £ <-k,iJa,i>/*(ri)/fo) • (10) 



1. 



Here, /(r^) is any normalized function J^ii l/( r 
This definition encompasses and generalizes the plane 
waves described in Refs. [16 , and allows us to calcu- 
late the entanglement witness W as a function of spatial 
separation, which may be of general interest outside this 
particular example. 

To probe spatial dependence, we calculate the entan- 
glement between two sets of lattice bins r s= i... m and 
r w=i...n separated by a distance Sr using the witness W 
with the function 



fin) 



1 if T{ G r 8 , 

exp(i0) if n G r w , 
otherwise. 



(ii) 



For given 5r, W can be minimized with respect to 0. 

Numerical results. — We illustrate this technique us- 
ing a ID chain of spin j = 1 atoms with n s = 200 sites 
and n a = 10 atoms per site [17 , which is related to 
the bilinear-biquadratic Hamiltonian, which has a rich 
phase diagram displaying ferromagnetic, critical, dimer- 
ized, and Haldane phases, each with distinctive spatial 
correlation signatures [T8H24] . We have also checked that 
our method can be used to prepare correlation signa- 
tures of more exotic quantum phases, such as the critical 
phase of the bilinear-biquadratic Hamiltonian, which has 
a structure factor peaked at k = ±2tt/3 [20], and alge- 
braically decaying correlations with characteristic period- 
3 oscillations [T9] . 

We demonstrate the preparation of spin correlations 
Cdes(5r) with: (a) an exponential decay exp(— r/£) with 
a correlation length £, corresponding to gapped phases, 
such as spin liquids which are conjectured to appear in 
the vicinity of high-T c superconductivity [2 , 25 ; and (b) 
an algebraic decay r - ^, corresponding to critical phases 
and quantum critical points of the phase diagram. We il- 
lustrate case (a) with £ = 5 and case (b) with £ = 0.7. We 
compute the f p corresponding to Cdes(Sr) as described 
above, apply the pulses in sequence to the atoms, and 
compute the resulting real-space spin-correlations 



coincide well with the calculated correlation function in 
both cases (first column). For case (a), the exponential 
decay is maintained over several orders of magnitude, 
and fits to the short-range behaviour yield a correlation 
length close to the desired £ = 5. For case (b), a clear 
algebraic decay is seen with a fitted ( w 0.4. Devia- 
tions from the desired parameters induced by finite opti- 
cal depth could be further compensated by adjusting the 
C p appropriately. The real-space correlation signature 
can also be extracted by fitting the Fourier transform 
of Cdes(^) to the covariances r(fc, — k) (second column), 
which are the observables that are measured and ma- 
nipulated in the experiment. Finally, we calculate the 
entanglement witness W for a single bin entangled with 
a chain of 106 bins (third column). In both cases, W 
is minimized for <j) = and decays exponentially (alge- 
braically) with Sr following the spatial behavior of the 
spin correlation function. Notably, the entanglement is 
stronger for algebraic decay. 

Outlook. — We have demonstrated that, with a sim- 
ple modification of the experimental scheme discussed 
in Ref. [4], it is possible to engineer a quantum lattice 
gas with an arbitrary non-classical spin correlation func- 
tion. We have illustrated the procedure with two ex- 
amples mimicking the quantum phases of the bilinear- 
biquadratic Hamiltonian, demonstrating how to prepare 
exponentially- and algebraically-decaying correlations. 
We have also checked that the method can be extended 
to the spatial correlations signatures of more exotic quan- 
tum phases. We have generalized the entanglement wit- 
ness proposed in Refs. [16^ and shown that the engineered 
spin-correlations entail multimode atomic entanglement. 
In our calculations, we make conservative assumptions 
about the experimental parameters, leaving considerable 
scope for further optimization of the procedure, which is 
readily extendible to higher dimensions and larger-spin 
systems, with both fermionic or bosonic atoms. 
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EE r -( f - r ^^ r )/ r(0) ( 12 ) 



i=l a 



after all pulses have been applied. The only remaining 
free parameter is then the maximum coupling strength 
max p {Cp}, which we set to « 0.95, ensuring that the ap- 
proximations suggesting the used coupling strengths ([9| 
are valid [26] . 

The numerical results are shown in Fig. [2] At large 
optical depth, the desired correlation signatures Cdes(5r) 



philipp.hauke@icfo.es 
[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008); M. Lewenstein, A. Sanpera, V. Ahufin- 
ger, B Damski, A. Sen, and U. Sen, |Adv. Phys. 56, 243] 
(2007) ; T. Esslinger, Annu. Rev. Cond. Mat. Phys. 1, 129 
(2010); M. Lewenstein, A. Sanpera, and V. Ahufinger, 
Ultracold Atoms in Optical Lattices: Simulating Quan- 
tum Many- Body Systems (Oxford University Press, Ox- 
ford, 2012). 

[2] P. W. Anderson, Science 235, 1196 (1987); G. Kotliar, 
|Phys. Rev. B 37, 3664 (1988)[ P. A. Lee, N. Nagaosa, 



5 



and X.-G. Wen, Rev. Mod . Phys. 78, 17 (2006)] 
[3] K. Eckert, L. Zawitkowski, A. Sanpera, M. Lewenstein, 

and E. S. Polzik, \Phys. Rev. Lett. 98, 100404 (2007)| 
[4] K. Eckert, O. Romero-Isart, M. Rodriguez, M. Lewen- 
stein, E. S. Polzik, and A. Sanpera, |Nat. Phys. 4, 50] 
(2008) 

[5] T. Roscilde, M. Rodriguez, K. Eckert, O. Romero-Isart, 
M. Le wenstein, E. Polzik, and A. Sanpera, |New J. Phys. 
11, 055041 (2009)| G. De Chiara, O. Romero-I sart, and 
A. Sanpera, |Phys. RevTT 83, 021604 ( 2011)] C. Weit- 
enberg, P. SchauB, T. Fukuhara, M. Cheneau, M. En- 
dres, I. Bloch and S. Kuhr, |Phys. Rev. Lett. 106, 
215301 (2011)1 G. De Chiar a and A. Sanpera, |J. Low] 
Temp. Phys. 165, 292 (2011) ; J. Meineke, J.-P. Bra ntut 
D. Stadler, T. Mull er, H. Moritz, and T. Esslinger, Nat] 
| Phys 8, 455 (2012)| 

[6J J. Appel, P. J. Windpassinger, D. Ob lak, U. B. Hoff 
N. K j aergaar d , and E. S. Polzik, [Proceedings of] 
the National Academy of Sciences 106, 10960 (2009) ; 
T. Takano, M. Fuyama, R. Namiki, and Y. Takahashi, 
|Phys. Rev. Lett. 102, 033601 ( 2009); M. H. Schleier- 
Smith, I. D. Leroux , and V. Vuletic, |Phys. Rev. Lett.] 
104, 073604 (2010)| Z. Chen, J. G. Bohnet, S. R. Sankar 
J. Dai, and J. K. Thompson, Phy s. Rev. Lett . 106, 
133601 (2011 ) ; R. J. Sewell, M. Koschorreck, M. Napoli- 
tano, B. Dubost, N. Behbood, a nd M. W. Mitchell, 
(2011), |a7Xiy:1111.6969 [quant-ph]! 

[7] G. Tot h and M. W. Mitchell, |New J. Phys. 12, 053007 
| (2010)| 

[8] K. Hammerer, K. M0lmer, E. S. Polzik, and J. I. Cirac, 
|Phys. Rev. A 70, 044 304 (2004); S. R. de Echaniz, 
M. Kosc horreck, M. Napolitano, M. K ubasik, and M. W. 
Mitchell, |Phys. Rev. A 77, 032316 (2008)| 
[9] M. Kosc horreck, M. Napolitano, B. Dubost, a nd M. W. 
Mitchell, |Phys. Rev. Lett. 105 , 093602 (2010)| 

[10] K. M0l mer and L. B. Madsen, |Phys. Rev. A 70, 052102 1 
(2004)1 L. B. M adsen and K. M0lmer7 |Phys. Rev. A 70] 
052324 (2004)] M. Koschorreck and M. W. Mitchell, J. 
Phys. B 42, 195502 (2009). 

[11] Alternatively, data with (J a (2k p )) < T (0) V could be 
post-selected based on the measurement outcomes. 

[12] The pulse with k = 2k p also introduces small correlations 
to the zero mode T zz (k, 0) ^ 0, which results in a change 
of F zz (k,k) and F zz (k,—k) by subsequent pulses with 
k' 7^ k. Below we demonstrate that this has negligible 



effect on the outcome. 

[13] This restriction could be relaxed with a more sophisti- 
cated strategy that takes into account the change in the 
correlations at one wavevector due to squeezing a differ- 
ent wave vector. 

[14] G. Gied ke and J. Ignacio Cirac, |Phys. Rev. A 66, 032316 

I (2002)| 

1 15 1 T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 
032110 (2002)1 A. Osterloh, L. A mico, G. Falci, and 
R. F azio, |Nature 416, 608 (2002)| R. Jozsa and N. Lin- 
den, [P roceedings o f the Ro yal Society of Londo n. Se-| 

Eries A: Mathematical, Physical and Engineering Sciences 
459, 2011 (2003)| F. Verstraete, M. Popp, and J. I. 
Cirac, Phys. Rev. Lett. 92, 027901 (200 4); O. Giihne 
and G. Toth, |Physics Reports 474, 1 (2009) 
[16] P. Krammer, H. Kampermann, D. BruB, R. A. Bertl- 

mann, L. C. Kwek, and C. Macchiavello, |Phys. Rev.| 

Lett. 103, 100502 (2009)| M. Cramer, M. B. Plenio, 
and H. Wunderlich, |Phys. Rev Lett. 106, 020401 (2011)| 
G. De Chia ra and A. Sanpera, | J. Low Temp. P hys. 165, 
292 (20TT)] 

[17] Note that the same results generalize to a single atom 
per site, as long as we bin the atoms into n s bins with n a 
atoms per bin, and redefine the coupling constant C p as 

an average over th e n a atoms in each bin. 

[18] F. D. M. Haldane, Phys Rev. Lett. 50, 1153 (1983)| 
[19] G. Fath and J. Solyom, |Phys. Rev. B 44, 11836 (1991)] 
[20] R- J- Bursill, T. Xiang, and G. A. Gehring, | Journal of| 
Physics A: Mathematical and General 28, 2109 (1995)| 
U. Schollwock, T. Jolicceur, and T. Garel, Phys. Rev. B 
I 53, 3304 (1996)| 

[21] A. ImambekovTM. Lukin, and E. Demler, [Phys. Rev. A 
I 68, 06 3602 (2003) 

|22| K. Buchta, G. Fath, O. L egeza, and J. Solyom, |Phys.| 

Rev. B 72, 054433 (2005)] 
[23] M. Rizzi, D. Rossini, G. De Chiara, S. Montangero, and 

R. Fazio, Phys. Rev. Lett. 95, 2 40404 (2 005)] 
[24] A. Lauchli, G Schmid, and S. Trebst, |Phys. Rev. B 74~ 
I 144426 (2006)| 

[25] G. Misguich and C. Lhuillier, "Frustrated spin systems," 
(World Scientific, Singapore, 200 4) p. 229; C. Lhuillier, 
(2005) , |arXiv:cond-m at /0502464] 

[26] This is a conservative choice: we could increase C p and 
calculate the f p taking into account the change in the 
covariances T zz (ki, fe) after each pulse in the sequence. 



